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Given an ensemble of mixed qubit states, it is possible to increase the purity of the constituent states using a 
procedure known as state purification. The reverse operation, which we refer to as dilution, produces a larger 
ensemble, while reducing the purity level of the systems. In this paper we find asymptotically optimal procedures 
for purification and dilution of an ensemble of independently and identically distributed mixed qubit states, for 
some given input and output purities and an asymptotic output rate. Our solution involves using the statistical 
tool of local asymptotic normality, which recasts the qubit problem in terms of attenuation and amplification of a 
single-mode displaced Gaussian state. Therefore, to obtain the qubit solutions, we must first solve the analogous 
problems in the Gaussian setup. We provide full solutions to all of the above, for the (global) trace-norm figure 
of merit. 
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I. INTRODUCTION 

When implementing any quantum information protocol, the 
states we wish to employ and manipulate are inevitably af- 
fected by decoherence effects, which diminish their purity 
and consequently their resource power. There exist several 
well-established methods to protect against such undesirable 
factors: strengthening the entanglement resource using dis- 
tillation methods [T| or employing a quantum eiTor correc- 
tion scheme |2| to encode our 'fragile' states into some larger, 
more unyielding system. The method we study in this paper is 
that of state purification fT '^], a procedure which takes as in- 
put an ensemble of identical copies of an arbitrary (unknown) 
state and produces as output a smaller ensemble of identical 
states with higher purity. This can be seen as a special case 
of the more general problem of inverting the effect of a noisy 
channel on ensembles of states, the channel being the depo- 
larising one in the present study. 

There already exists several theoretical results for purifi- 
cation of n i.i.d. (independently and identically distributed) 
mixed qubits, notably Refs. |l3]|4l, where optimal purification 
algorithms for various formulations of the purification prob- 
lem are provided. Purification of an ensemble of mixed qubit 
states has also been found to occur in the context of 'super- 
broadcasting' lis), an 71 — ^ m cloning procedure which can 
actually result in purified clones for n > 4 and sufficiently 
mixed input states (the noise present is merely shifted from lo- 
cal states into correlations between output states). For n > m, 
superbroadcasting is actually equivalent to the optimal purifi- 
cation procedure of Ol. Experimentally, purification has been 
achieved in |6|, which implemented the methodology of ||3j 
and demonstrated optimal purification for the case n = 2. 

Beyond the entanglementology (phenomenology of entan- 
glement), judging the performance of a purification protocol 
requires a figure of merit (FoM) which measures the departure 
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from the ideal transformation. Two types of FoM have been 
considered in the literature, with very different results. The lo- 
cal FoM is built upon the comparison of the reduced states of 
individual output systems with the target state. In this case, a 
complete reversal of the depolarising channel may be obtained 
asymptotically with the size of the input ensemble, and with 
arbitrarily high output rate m/n H. The global FoM com- 
pares the joint state of the output with that of a product of inde- 
pendent target states. This is a more demanding criterion. For 
example if the output systems are independent and identically 
prepared then the global fidelity scales as m) = i^™ 
where < 1 is the fidelity of an individual output state with 
respect to the target state. Indeed, it has been shown I?) that no 
protocol can achieve asymptotic purification [F{n, m) — ^ l) 
to pure target states at a finite rate m/n. The global figure or 
merit is relevant whenever we deal with the collective state of 
the output rather than the individual constituents, as in the case 
of state transfer between atomic ensembles and light. Addi- 
tionally, it can serve as a "measure of coiTelations" when the 
individual constituents of the output states are known to be 
exactly in the target state, as in superbroadcasting. This hy- 
pothesis will however not be pursued in this paper 

The above no-go theorem motivates us to consider the ques- 
tion whether the depolarising channel can be reversed with a 
positive asymptotic output rate, when the target states (i.e. the 
states prior to applying the depolarising channel) are mixed. 
We show that this is indeed possible, and compute the maxi- 
mal purification rate for given input and target purity, and the 
optimal FoM for approximate purification at a fixed rate which 
is higher than the maximal one. 

We also consider the opposite process of dilution in which, 
starting from an ensemble of n identically prepared states, 
we produce a larger ensemble consisting of m independent, 
but more mixed states. Dilution shares similarities with the 
process of optimal n — > m quantum cloning Q, but while 
in cloning the rate m/n is fixed, and one aims at generating 
clones as close as possible to the input states (with respect 
to a local or a global FoM), in a dilution procedure we set a 
target level of output purity and look for the optimal rate for 
generating such target states. 
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qubit problem 


Gaussian problem 


state model 


ensemble of n i.i.d mixed qubits 
[Eq. {T6|] 
71 2> 1: number of copies 
r with ||r|| < 1: Bloch vector; 


single-mode displaced Gaussian state 
[Eq. 0] 
Q £ C: displacement; 
s G (0, 1): purity parameter 


input 






target 


^ XrQ+ku/ ^/m 


Kl 


procedure 


purification 

A > 1, m < n 


attenuation 

S2 < Si, fc < 1 


procedure 


dilution 

A < 1, m > n 


amplification 

S2 > Si, k > 1 



TABLE I: Summary of the notation adopted in the present paper. For the qubit problem, we aim at optimising the output-vs-input rate 
m/n (maximising it for purification, and minimising it for dilution) at given input Bloch vector ro + u/^/n and scale factor A. For the 
corresponding Gaussian problem, we aim at finding the maximal value of the displacement ratio k, such that attenuation or amplification can 
be realised perfectly, for given target temperature parameters Si and S2 and unknown displacement a. The framework of local asymptotic 
normality provides a rigorous link between the two problems, as for n ^ 1 the local Bloch vector u is mapped into the displacement a of a 
single-mode coherent thermal state. 



In deriving the asymptotic results, the key mathematical 
tool is that of local asymptotic normality (LAN), a fundamen- 
tal 'classical' statistics technique |8| which was recently ex- 
tended to the context of quantum statistical models ||9HT2). 
In the quantum case, LAN dictates that the collective state 
of n i.i.d. quantum systems, can be approximated by a joint 
Gaussian state of a classical and a quantum continuous vari- 
able (CV) systems. This has been used to derive asymptot- 
ically optimal state estimation strategies for mixed states of 
arbitrary finite dimension and also in finding quantum 
teleportation benchmarks |14| and optimal quantum learning 
procedures 1 15 1 for multiple qubit states. The general strat- 
egy is to recast statistical problems involving n i.i.d. quantum 
systems into the simpler setting of Gaussian states. The op- 
timal solution for the corresponding Gaussian problems can 
then be used to construct asymptotically optimal procedures 



for the original one. In section III we sketch how this could 



be physically implemented, and more details can be found in 

Eol. 

Following this methodology, we transform the qubit purifi- 
cation and dilation problems into those of optimal attenuation 
and amplification for a one-mode CV system in a Gaussian 
state, together with a classical real-valued Gaussian variable, 
both with known variance but unknown means. In attenuation 
we reduce the variance of a displaced Gaussian state, at the 
price of simultaneously reducing its amplitude, while in am- 
plification we increase the amplitude, as well as the variance. 
For both problems we use a FoM based on maximum trace- 
norm distance, and show that the optimal attenuation chan- 
nel is obtained by applying a beamsplitter, while the optimal 
amplification is implemented by a non-degenerate parametric 
amplifier A similar scheme for the attenuation of Gaussian 
CV states has been proposed and experimentally implemented 
in lfT6i . Parametric amplification has been investigated in 1 17 - 
[T9l , and demonstrated experimentally in |20|. In particular, 
the same amplifier is optimal for a FoM based on the mini- 



mum amount of added noise ifTTlfTSl. However, whilst these 
transformations are well known candidates for our protocols, 
to the best of our knowledge a proof of their optimality with 
respect to the FoM chosen in this paper had not been obtained 
in the literature. Our proof relies on a covariant channels op- 
timisation technique developed in fT4"21|. We find that for 
given input and output purity parameters, there exists a range 
of values for the ratio k between output and input displace- 
ment, such that attenuation or amplification can be realised 
perfectly, and we compute the maximal (optimal) value k^, 
as a function the two purities. In the parameter range where 
the procedures cannot be accomplished perfectly, we give the 
exact expression for the optimal FoM. 

A schematic summary of the problems addressed in this pa- 
per is provided in Table |l] The paper is organised as follows. 
In Section [lljwe formulate and solve the two quantum Gaus- 
sian problems, and the corresponding classical one. In Sec- 
tion III we use this result in conjunction with LAN to find 



asymptotically optimal purification and amplification chan- 
nels for states of n i.i.d. mixed qubits. We draw our con- 



cluding remarks in Section IV The proofs are collected in 
Appendix A. 



II. OPTIMAL ATTENUATION AND AMPLIFICATION OF 
GAUSSIAN STATES 

A. Classical Case 

Before we move onto the quantum case, it is instructive and 
relevant to consider the corresponding problems for classical 
random variables. In the classical scenario, the analogue of 
'attenuation' ('amplification') is a procedure which reduces 
(increases) the mean and variance of a given random variable. 
The analogue to our quantum problem would then be to find 
a transformation K which maps a real-valued normally dis- 
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tributed random variable X ^ N{u, Vi) of arbitrary mean u 
and fixed variance Vi, into a variable Y ^ A^(fcw, V2) such 
that the risk 

R„,UK;Vi,V2,k) ^2svip\\K{Niu,Vi)) - N{ku,V2)U 

U 

(1) 

is minimised. Here k represents a fixed constant, where 
< k < 1 means attenuation and A: > 1 means amplification 
of the Gaussian variable X, and we choose the interesting case 
where Vi > V2 in the case of attenuation, and Vi < V2 for am- 
plification. The notation ||P - Q||tv = sup{|P(A) - Q{A)\ : 
A G J^}, for the cr-algebra represents the total variation 
distance between the probability distributions P and Q which 
reduces to one-half of the Li-distance between their proba- 
bility densities in the case of mutually absolutely continuous 
distributions |22|. 

The solutions of both classical and quantum versions of 
this problem rely on the notion of 'covariance'. Consider the 
transformation 



X ^ K{X) = kX 



(2) 



where X and Z are independent random variables, Z having a 
fixed variance and vanishing mean. Such a (classical) channel 
is covariant, in the sense that 



K{X + C)= K{X) + kC 



(3) 



for any constant C. Such transformations can be shown to 
not only minimise ([T]l, but also to render it independent of 
expectation so that the FoM becomes 

i?max(i^; Vu V2, k) - 2|| A'(iV(0, VtS) - iV(0, F2)|ltv. 

It is easy to see that if 



k<k\^\Vi,V2) := 



(4) 



then the target distribution can be achieved exactly, with the 
appropriate amount of Gaussian noise in the variable Z. As 
we shall see in the next section, there exists an analogous 
range (12 1 for the quantum Gaussian transformation. 

As for the case k > k^'^^Vi, V2), it can be shown 12211 that 
the optimal choice for Z in (|2]) is Z = 0, as one would expect, 
so that the optimal figure of merit is 



n 

K 



1 



,{K-VuV2,k) 



yW2 



dx. (5) 



Henceforth, we will denote by K* the optimal transforma- 
tion for the two cases discussed above. 



B. Quantum Case 

In this Section we consider the following: given a Gaussian 
state $0, of a one-mode CV quantum system, with known co- 
variance and unknown displacement a, we would like to op- 
timally attenuate (amplify) it, that is transform it into a state 




V 



(a) 



(b) 



FIG. 1 : (Color online) Schematic phase-space diagram for (a) atten- 
uation and (b) amplification of a displaced Gaussian state <l>^i . 



with smaller (greater) covariance and displacement ka, with 
the largest possible proportionality constant k. Let 

Wa '■— exp(aa^ — aa) 

denote the Weyl operators where a e C and a, the creation 
and annihilation operators satisfying [a, a^ — 1 and 



a\n) 



I) 



n > 0, 



where {|n)},i>o is the Fock basis of the Hilbert space H. For 
< s < 1 we denote by the centred, phase invariant 
Gaussian state 



n=0 



^){r 



(6) 



and by displacing it we obtain the family of Gaussian states 

Wc^^'Wl (7) 

Given two different mixing parameters si > S2 (si < S2) and 
a positive parameter A: < 1 (fc > 1) we would like to find the 
optimal attenuation (amplification) channel which maps the 
state close to the state for an arbitrary displacement 
a (see Fig. [T]). For any channel P : Ti{'H) we 
define the FoM called the maximum risk 



-Rmax(-P;Sl,S2,A:) 



and the minimax risk 



sup 



^? I 



R 



mm max 



(si,S2,fc) 



inf i?inax(-P;si,S2,A:). 



(8) 



(9) 



A channel is called 'minimax' if its maximum risk is equal to 
the minimax risk. We will show that (up to a trivial adjust- 
ment for a certain range of fc's) the optimal solutions to the 
attenuation and amplification problems are, respectively, the 
beamsplitter and parametric amplifier. 

We start by defining a specific channel denoted in both 
cases P*, then show that it is optimal and compute the mini- 
max risk. For si > S2 and fc < 1, the attenuation channel is 
implemented by the action of a beamsplitter with reflectivity 
k acting on an input mode a prepared in a state $^1, and a 
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FIG. 2: (Color online) FoM's for optimal attenuation (left) and optimal amplification (right) of displaced Gaussian states. Top row: (a) Plot 
of the minimax risk -Rmax [Eq. versus k for the optimal attenuation procedure P*, where si = 0.8 and S2 = 0.4; the optimal value kg 
is highlighted with a big dot on the graph, (b) Plot of the minimax risk -Rmax [Eq. l |14^ ] versus k for the optimal amplification procedure P* , 
where si = 0.4 and S2 = 0.8; the optimal value ko is highlighted with a big dot on the graph. Bottom row: 3D Plots of the minimax risk 
i?max versus the parameter k and the output/input temperature ratio S2/S1 with si = 0.5, for (c) the optimal attenuation procedure, and (d) 
the optimal amplification procedure; in both plots, the thick curve depicts kg as a function of S2/S1. 



second ancillary mode b in the vacuum state uj ~ |0) (0|. The 
output mode c of the channel is 

c = k^a+ y/l- k^b. (10) 

For S2 > si and fc > 1, the channel is a parametric amplifier, 
whose action is represented by the following transformation 
on the input mode a and an ancillary mode b prepared in the 
vacuum state: 

c = ka+^/k^ - 16^ (11) 

We note that for each pair (si, S2) there exists a range of 
parameters k for which i?minmax(si, S2, k) = 0, i.e., the pro- 
cedures can be accomplished perfectly. Indeed it can be easily 
verified that, for k given by 

;att/ N /s2(l-Sl) .amp, n /l - 

(12) 

the channels ( [TOj l and respectively ( [TT| ) produce exactly the 
target state "I"^^. Moreover, if /c < /cq then the output of P* 
is the state with s < S2, and the target can be still per- 
fectly achieved by adding an appropriate amount of Gaussian 
noise. For later use, when fc < fcp we will denote by the same 
symbol P* this modified channel. From now on we consider 



the less trivial situation fc > fco, corresponding to the regime 
where perfect amplification or attenuation are impossible. We 
then state the following theorem and lemma, whose proofs are 
given in Appendix A: 

Theorem II.l. If k < k^ then the minimax risk for attenua- 
tion (amplification} is zero. Ifk> fco, the minimax procedure 
is P*, i.e. the beamsplitt er (|7Q| ) in the case of attenuation, or 
the parametric amplifier ( |i7| ) in the case of amplification: 

(si,S2,fc). (13) 

Lemma II.2. Ifk <kQ, then the minimax risk for attenuation 
(amplification) is given by 

i?mi„max(si, S2, fc) = 2(5"» + l - S^" + ^), (14) 

where mo is the integer part of 

ln[(l-S)/(l-S2)]/ln(s2/5), 
and s takes the values 

.Sifc^ 1 — Si 

= l-,si+.ifc2' ^""^ = 1 - ^^^^ 
in the case of attenuation and respectively amplification. 

The risk for both processes is plotted in Fig.[2][(a)-(d)]. In 
Figure [3] we plot fcp for both processes as a function of the 
input and output purity parameters si and S2. 



5 



■^2 



1.0 
0.8 
0.6 
0.4 
0.2 



0.0 



1.0. 

0.8 I 

0.6 
0.4 
0.2 
0.0 



0.0 0.2 0,4 0.6 0.8 1.0 

(a) 




0.0 0.2 0.4 0.6 0.8 1.0 
(b) 



single system. Note also that the risk at a fixed point r can 
always be made equal to zero by simply preparing the target 
state for that point. To take into account the overall perfor- 
mance of a procedure, one can either integrate the risk with 
respect to a prior distribution over states (Bayesian statistics) 
or take the maximum over all states (frequentist statistics). We 
adopt the latter viewpoint, and in addition we will consider a 
more refined version of the maximum risk called local maxi- 
mum risk around Tq 



RmaxiQn', Tq, A) 



sup i?((3„;r, A). 



(18) 



FIG. 3: (Color online) Contour plots of (a) ko versus the purity pa- 
rameters si and S2 for the optimal attenuation procedure, and (b) 
ko~^ versus the purity parameters si and S2 for the optimal amplifi- 
cation procedure. 



III. ASYMPTOTICALLY OPTIMAL PURIFICATION AND 
DILUTION FOR ENSEMBLES OF QUBITS 

We turn now to the problem of finding optimal purification 
and dilution schemes for ensembles of identical qubits. We 
denote by the qubit state with Bloch vector r = {rx,ry,rz) 



1 . 



(16) 



where rcr = r^a^+ryay+rzO-z and (7^ are the Pauli matrices. 
We are given n identical qubits prepared in the state pr and we 
would like to produce m identical qubits in the state p\r with 
m as large as possible, for a fixed positive parameter A. When 
A > 1, the aim is to "purify" the state, and when < A < 1 
we want to "dilute" the state with the benefit of obtaining more 
copies. Clearly, for purification the output state is physical 
only if A satisfies A||r|| < 1. This can be achieved by letting 
A depend on r, or by restricting to those input states which 
satisfy the property. To illustrate the latter, suppose we would 
like to reverse the action of the depolarising channel 

C : Pr ^-> Pr/A, 

then the input states of the purification channel automatically 
satisfy the requirement. As to the former, our asymptotic anal- 
ysis will produce a local FoM which only depends on the 
value of A at a particular state, so for simplicity we will as- 
sume it to be constant. 

A purification (dilution) procedure is a quantum channel 



M{C^ ) 



mapping n-qubit states to m-qubit states, and its performance 
is measured by the FoM (risk) 



i?(Q„;r,A) := HQnipf^ - pTh 



(17) 



Note that this is a global rather than a local risk, in the sense 
that it measures the distance between the output and the joint 
product state, instead of comparing their restrictions to each 



In asymptotic statistics the local maximum risk is more infor- 
mative that the 'global' one since it captures the behaviour of 
the procedure around any point in the parameter space, rather 
than that of the worst case. The radius of the ball over which 
we maximise is slightly larger than the precision of 
with which we can estimate the state parameters, so that the 
definition of the local risk does not amount to assuming any 
prior information about the parameter Indeed one can use a 
small sample n^^^ <C n of the input systems to obtain a rough 
estimate of the Bloch vector r such that the obtained estimator 
To will be in a ball of size n^^/^^'^ around r, with probability 
converging to one as n cx). With this additional informa- 
tion, one can then apply the purification (dilution) channel to 
the remaining systems, with no loss in the asymptotic opti- 
mal risk (see below). The local maximum risk is a standard 
FoM in asymptotic statistics and it is has been used in quan- 
tum statistics in |9 , 10 15 1 to which we refer for more details, 
and for its relation to Bayesian methods. 

Up to this point the number of input and output systems n 
and m were fixed, with n considered to be large. However, for 
a non-trivial asymptotic analysis, m should be an increasing 
function of n, more precisely we consider the optimal purifi- 
cation (dilution) procedure for a fixed rate 

m(n) 

A = lim — ^ > 0. 

n— f oo n 

Indeed from our fixed rate analysis it can easily be deduced 
that in the case of a sub-linear dependence m{n) = o{n), one 
can produce m output copies of arbitrary purity with vanish- 
ing local maximum risk. On the other hand, by similar reason- 
ings, one may expect that if m{n)/n is unbounded, then the 
best strategy should be to estimate the state and reprepare m 
independent copies of the estimator ('measure and prepare' 
strategy 1231 ). We leave this statement as a conjecture, and 
from now on we will assume that the rate A is given and fixed. 
For any sequence Q :— {Qn} of procedures we define the 
asymptotic local maximum risk at Tq by 

i?(Q;ro,A,A) := limsupi?max((3n;ro, A), (19) 

and we would like to find an optimal (minimax) strategy 
whose asymptotic risk is equal to the minimax risk 

-Rminmax(ro, A, A) : = lim sup inf i?i„ax(Q« ; Tq , A) . (20) 

n— ^OO Qn 
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In other words, we will answer the following question: for 
given purification (dilution) constant scale factor A and input- 
output rate A, what is the minimax risk i?mmmax(ro, A, A) 
and which is the procedure that achieves it? In particular, we 
will find that the minimax risk is zero for a range of parame- 
ters (ro, A, A), and we will identify the maximum value Aq""^ 
(Aq'') for which the purification (dilution) can be performed 
with asymptotically vanishing risk. These rates are the qubit 
analogues of the constants fco defined in ( [T2] i. 

The main technical tool is the theory of local asymptotic 
normality (LAN) developed in |9 - 12| as an extension of a key 
concept from (classical) asymptotic statistics llsl [24l . LAN 
means that the joint quantum state of identically prepared 
(finite-dimensional) systems can be approximated in a strong 
sense by a quantum-classical Gaussian state of fixed variance, 
whose mean encodes the information about the parameters 
of the original state. In this way, a number of asymptotic 
problems can be reformulated in terms of Gaussian states, for 
which the explicit solution can be found, e.g. state estimation 
II25I . teleportation benchmarks lfT4l . quantum learning ifTSl . 
system identification |26|. For the purposes of this paper we 
give a brief description of LAN for mixed qubit states. Let 



Pro+a/^/n 



1 



1 + (ro + u/Vn)cr) 



denote a qubit state in a the neighbourhood of a fixed and 
known state prg, which is uniquely characterised by an un- 
known local parameter u. The family of n-qubit states 



(21) 



will be called the local statistical model at Tq. Additionally, 
we define a classical-quantum Gaussian model 

AA:= {Ar„«)$„ :ueM^} (22) 

where N„ := N{uz, 1 — ||ro|p) is a normal (Gaussian) distri- 
bution on M with mean and variance 1 — ||ro|p, and 



|ro| 



2||rn 



1 



|rol 



is a displaced thermal Gaussian state of a one-mode CV sys- 
tem (cf Section |IIB[ ) with known covariance matrix charac- 
terised by the purity parameter s (with zero squeezing) and 
unknown means proportional to {ux,Uy). Now, the mathe- 
matical statement of LAN H is that there exist two sequences 
of channels T = {r„} and S = {5'„} with 

T„ ; M{C^") ^ L\R)(g>Ti 
Sn : L\R)®Ti -> Af(C2") 



such that 



lim sup \\T.n{p:)-N,^<i>J 
lim sup ||p;'-5„(Aru®$u)|| 

||ui|<n- 



: 0. 



In the above formulas, 7i is associated to the trace-class op- 
erators of the CV system, and the norm-one || • ||i denotes 
respectively the trace-norm for the quantum part and the Li- 
norm for the classical part. The physical implementation for 
the channels T„ and Sn, detailed in 1 13 1, is realised via a spon- 
taneous emission coupling of the n qubits to a Bosonic field, 
and subsequently letting the qubits Teak' into this environ- 
ment. Since there is no correlation between atoms and field, 
the statistical model decouples into a Gaussian state $„ asso- 
ciated to the field, and a classical statistical mixture of atoms, 
distributed according to N^. The corresponding operations of 
attenuation and amplification may then be carried out in the 
field in the optimal way. 

In our case we need to consider mixed qubit states, which 
means that the collective state has non-zero components in 
all irreducible representations of SU(2) (all values of the to- 
tal spin). In fact the traces of the different blocks of given 
total spin form a probability distribution which (after centring 
and scaling) converges to the classical Gaussian component of 
the limit model in LAN. A typical block state of definite to- 
tal spin can be mapped isometrically into the Fock space of a 
one mode CV system, and converges to the quantum Gaussian 
component of the limit model. This transfer can be imple- 
mented in principle by a creation-annihilation coupling with 
a Bosonic field in which the state 'leaks" after a short time. 
The classical part (total spin) can be measured by coupling 
subsequently with another Bosonic field, and performing an 
indirect measurement of in the field. Since after the first 
step, all blocks are in the state, a measurement of is 
implicitly a measurement of the total spin. 

The above convergence can be interpreted as follows: the 
quantum statistical models Vn can be mapped into the Gaus- 
sian model Af and vice-versa, by means of physical operations 
(quantum channels) with vanishing norm-one error. From the 
statistical point of view, in many situations this convergence 
is strong enough to allow us to map a statistical problem con- 
cerning the model Vn to a similar one concerning the simpler 
model Af. 

In the case of purification or dilution of qubits, the map- 
ping into a Gaussian problem is illustrated in the diagram be- 
low. We first give a detailed description of the steps involved, 
and then prove that our procedure is optimal (asymptotically 
minimax). 



^'0 + :^ 



p 



Aro + 

Sm 



(23) 



K* 



Step 1: Localisation. We are given n identical qubits in 
an arbitrary mixed state p^. We measure a small proportion 
n^^'^ <C n of the qubits, to obtain a rough estimator of 
the state. By standard concentration results, with asymptoti- 
cally vanishing probability of error, the actual state is in a local 
neighbourhood of of size ?t,~^/^+'^, so that the remaining 
qubits can be parametrised as in the local model ( |2T] i. In the 
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FIG. 4: (Color online) Schematic Bloch-sphere geometry for (a) pu- 
rification and (b) dilution of qubits. A change in Bloch vector magni- 
tude by factor A is reflected by a change in the corresponding Gaus- 
sian state displacement. The two are then related by A = k^/njrn. 



same time, the target single-system output state p\r belongs 
to the local model 



Aro+- 



with local parameter (see Figure [4]i 



u' = AW-u = AAi/2u:=fcu. 



Step 2: Transfer to the Gaussian state. We apply the map 
T„ to the qubits and obtain a classical random variable and 
a single-mode CV system whose states are approximately 
Gaussian [see (|22])] 



|ro| 



Wo^^'Wl (24) 



Step 3: Optimal Gaussian purification (amplification). 
Since Q„j is a local model around Aro, the corresponding pa- 
rameter of the associated Gaussian model is 



S2 := 



l-A||ro| 
l + A||ro| 



and the family of Gaussian states is 

N{ku,, 1 - llfcrolP) ® Wkc^^'^Wl^. (25) 



In this step we attenuate (amplify) the Gaussian state p4| i in 
order to map it into, or close to ( p5] l, as described in Sectionjll] 
This means that we apply the optimal channel K* defined 
in section II A to the classical component N-^, and the opti- 



mal quantum attenuation (amplification) channel P* defined 



in Theorem II. 1 to the quantum Gaussian component 

Step 4: Mapping back to the qubits. We apply the channel 
S.„i to the classical variable and the output of the attenuation 
(amplification) channel to obtain a state of m qubits in the 
neighbourhood of the state pAro ■ 

By composing the channels employed in steps 2-A we ob- 
tain the overall channel 



Sn o [K* (g, p*) o r„ 





\ *o^"' > ko 


ko^'^ < ko 





0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 5: (Color online) Parameter range defining the different regimes 
of qubit dilution. In the shaded region below the boundary curve 
A = A(||ro||) [Eq. ji?}], the threshold k^g^ for a nonzero classical 
Gaussian risk is smaller than the corresponding threshold fc™^ for a 
nonzero quantum Gaussian risk, and Cases 3 and 4 of Theorem |III.l| 
apply for determining the optimal FoM for the qubit problem. Above 
the boundary curve, the situation is reversed, and Cases 2 and 4 of 
Theorem |III. 1 [ apply instead. 



which will be shown to be optimal. Recall that for fc < fco 
[see Eq. ( 12 1], the quantum component of the Gaussian 
target state can be prepared exactly. The same is true for the 



(c) 

classical component when k < , where 



' l-A^||ro| 
l-||rolP 



(26) 



is obtained by substituting Vi 
for the variances in (Hh. 



l-||ro||2,F2 = l-A2||ro|| 



This means that the total risk has different expressions over 
the following three intervals: it is zero when < fc < 
min(/co, ^0°^), it has one classical or quantum contribution for 

min(A:o, feQ^"*) < k < max(fco, ^o'^') and has both quantum 

and classical contributions for k > max(/co, ^o'^'')- For purifi- 
cation (corresponding to Gaussian attenuation), the ordering 
< /cq'* < fc'"^^ < 1 always holds, so the middle interval has 
a quantum contribution. However, for dilution (correspond- 

(c) 

ing to Gaussian amplification), the ordering of /cq and k^ ' 
depends on the parameters ||ro|| and A. In particular, we see 
the appearance of a boundary which demarcates the two sep- 
arate regimes of dilution, each defined by whether classical or 
quantum contributions to the risk take place first (see Fig.Bll. 
Namely, k^ 



< k^"'P for 



A<A(||ro||) = min 1 



i-lkol 

llroll 



(27) 



and kjf' > /cq'"'^ otherwise. 



Notice that inequality ( 27 1 is 



always satisfied for ||ro|| < 1/2 for all values of A. 

The relation between the output qubit rate A = m/n and 
the constants A, k can be inferred from the geometry of the 
Bloch sphere (see Fig.|4|) 



A(||ro||,A,fc) 



fc2 
A2' 



(28) 
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FIG. 6: (Color online) Optimal input (n) vs output (m) rates Ao 
of qubit production for the processes of perfect purification and dilu- 
tion, [Eq. 1I29I] , plotted versus the Bloch vector lengths before ( 1 1 Tq 1 1 ) 
and after (A||ro ||) the protocols. The left (right) side of the three- 
dimensional surface, corresponding to the region A > 1 (A < 1), 
represents the optimal rate m/n (n/m) for mixed qubit purification 
(dilution). 



In particular, the maximum output rates for which the asymp- 
totic risk is zero are obtained by setting: k = fcg** for purifi- 



cation, and 



for dilution. 
Explicitly, 



A< A(||ro||). 



ArdlrolU) 



fcn™P, Otherwise, 



A-^-||ro| 



A2(l 



(29) 



Af(||ro||,A) = 



||ro||+l/A 
•^^(Ikoll + l)' 



A< A(llroll). 



otherwise. 



The optimal rates ( p9| ) are plotted in Fig. [6] 

We can now state main result of this section whose proof is 
given in Appendix A. 

Theorem III.l. The sequence of purification (dilution) maps 

Ql o [K* ® P*) o r„ (30) 

is locally asymptotically minimax. 

We distinguish four cases for the minimax risk. 

Case 1." Zero risk. If k < min(fco, k^^^) then the minimax 
risk is zero. The optimal output rate Ag is given in \29) . 

Case 2: Quantum contribution. If k^ < k < k^ , then the 
purification ( dilution) minimax risk at rg is equal to the risk of 
the optimal Gaussian attenuation (amplification) scheme (|i4|) 



R 



ininmax 



(ro,A, A) = R, 



mmmax 



{k, 81,82). 



(31) 



Case 3: Classical contribution. If k\^^ k < kg, then the 
dilution minimax risk at rg is equal to the risk of the optimal 



classical Gaussian amplification scheme (|5|.- 



-^minmax (''0 7 A, A) -Rminmax 

iVi,V2,k). (32) 



Case 4: Classical and quantum contributions. If k > 
niax(/coj ^o"^'), then the purification (dilution) minimax risk 
is 



R 



minmax 



(ro,A,A) 
exp 



2fe2(i_||roP) 



(1 - S)F 



exp 



27rfc2(l - ||ro||2 



27r(l-A2||ro| 



(33) 



where s takes the values given in \15\ in the case of attenua- 
tion (for qubit purification) and amplification (for qubit dilu- 
tion) respectively. 

The optimal minimax risk for purification and dilution of 
qubits is plotted in Fig.|7]as a function of k. 



IV. CONCLUSIONS 

We have solved the practically relevant problem of optimal 
attenuation and amplification of displaced Gaussian states, 
with respect to the maximum norm-one distance FoM. As ex- 
pected, the optimal channels are implemented by the beam- 
splitter and parametric amplifier respectively, where the an- 
cillary state is provided by the vacuum in both cases. This 
solution was then used in conjunction with LAN, to construct 
optimal purification and dilution channels for ensembles of 
i.i.d. qubits as formulated in Theorem |III.l[ 

In the Gaussian case, we give an explicit expression of the 
FoM as a function of the variance parameters si and 82 of in- 
put and output states and the attenuation (amplification) factor 
k. In particular we identify the optimal value fco(si, 82) for 
which the protocol achieves the target state exactly. Similarly, 
in the multiple qubits case, we derive the FoM as a function 
of the input and output purity and the asymptotic input/output 
rate A, and identify the optimal rate Aq for which the the pro- 
tocol achieves the target collective state exactly. Both classi- 
cal and quantum Gaussian contributions to the risk have to be 
taken into account to calculate the maximum rates, and to pro- 
vide the optimal FoM for purification and dilution of multiple 
qubits, in the parameter range where the procedures cannot be 
accomplished perfectly. 

An interesting future project is to extend the techniques 
used in this paper to tackle the general problem of asymptot- 
ically optimal channel inversion for arbitrary channels on fi- 
nite dimensional systems. Such a study may provide efficient 
strategies to counteract the effect of various types of noise and 
decoherence processes, beyond the depolarising channel con- 
sidered in the present work. 
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0.0 0.2 0.4 0.6 O.i 
k 

(a) 




0.6 
0.5 
0.4 
0.3 



B 

0.2 



0.1 
0.0 



1.0 1.2 1.4 1.6 1.8 2.0 
k 

(c) 



FIG. 7: (Color online) Plot of the minimax risk -Rmax(Q*) [Eqs. i31 33|] for (a) optimal purification with ||ro|| — 1/3, A||ro|| — 4/5, (b) 
optimal dilution with ||ro|| = 4/5, A||ro|| — 1/3, and (c) optimal dilution with ||ro|| — 1/2, A||ro|| = 1/8, of n qubits, as a function of the 
local scale factor k. In panels (a) and (b), the risk is zer o for k < ko (black line), it is given by the quantum Gaussian FoM for ko < k < k^'^^ 
(blue line, corresponding to Case 2 of Theorem III. 1 



and it is then given by the joint quantum and classical contributions for k > k^ 



(red line, corresponding to Case 4 of Theorem |II .1^ . In panel (c), which describes a dilution characterised by a choice of parameters within 
the shaded region of Fig. dI we have instead that the risk is zero for k < fc*"' (black line), it is given by the classical Gaussian FoM for 
kf)^^ < k < ko (green line, corresponding to Case 3 of Theorem III.l >, and it is then given by the joint quantum and classical contributions for 
k > ko (red line, corresponding to Case 4 of Theorem [niTl. 
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Appendix A: Proofs 

Proof of TheoremllLll 

Proof. As the proof follows the lines of similar results in 
lfT4ll2n we will briefly sketch the main ideas. A covariance 
argument [27 , 28 1 shows that one may restrict the attention to 
channels which are displacement-covariant, in the sense that 
P(Wj$W|) = W^P($)W| for any input state p. For such 
channels the risk is independent of a and 

In the case of attenuation {k < 1) such channels are described 
by the linear transformation 

Catt = ka + \J\ — k'^b 

where a is the input mode, Catt is the output and b is an an- 
cillary mode prepared in a state r. Since the channel is com- 
pletely characterised by the state r, we will denote it by Pr- 
Similarly, for amplification (fc > 1) the output of the channel 
Pt is the mode 

Camp = ka+ \Jk'^ - 16^ , 



with h prepared in the state t. By a second covariance argu- 
ment with respect to phase rotations, and taking into account 
that $ is invariant under phase rotations, we obtain that r can 
be taken to be phase-invariant, i.e. it is a mixture of Fock 
states T = Ti\i){i\. In this case the output state will be di- 
agonal in the Fock basis and we write P^((f>*i) = 
and in particular corresponds to the output state when the 
ancilla is the vacuum. Similarly, we denote the coefficients 
of the Gaussian state <I>*i and by -pi — (\ — s\)s\ and 
qi = (1 — S2)s2- The proof reduces now to showing that, for 
any r, 

Ib'-glli > IIp^-^IIi- (ai) 

The key to proving this statement is the concept of stochastic 
ordering, whose definition we recall: 

Definition A.l (Stochastic Ordering). Let p — {pi : I e N} 

and q = {qi : I € N} be two probability distributions over N. 
We say that p is stochastically smaller than q {p di l) if 

m m 

^Pi>5l9i, Vm>0 (A2) 

1=0 1=0 

The following lemma holds for both the purification and the 
amplification scenarios: 

Lemma A.2. For any state t the following stochastic order- 
ing holds 

p'^ dp"- (A3) 

Proof. We treat the attenuation and amplification separately 
but the idea is the same in both cases: we reduce the state- 
ment about stochastic ordering to a simpler one where the in- 
put mode is in the vacuum. 
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Attenuation. We write the input mode as a = cosh(t)ai + 
sinh(t)a2 with ai,2 two fictitious modes in the vacuum state, 
and tanh^(i) :— si, which ensures that the state of a is 
Let t be such that sinh(t) k sinh(i) and denote 



The output mode is now 

Camp = smh{i){Rb^ + Tal) + cosh(f)ai 
= sinh(f)6^ + cosh(f)a 



T : 



(l-P) 



cosh2(t) 



R 



cosh(i) 



Then Cpur ^ ka + Vl — k'^b can be written as 



2 



cosh(f)(Tai + Rb) + sinh(?)a 
cosh(i)& + sinh(t) 



where b := Tai + Rb and 02 was relabelled d. The state of 
the mode b is given by 



00 
p=0 



\p){p\ 



so b is in the vacuum state if and only if b is in the vacuum. 
Thus it suffices to prove the stochastic ordering statement for 
the mode Catt written as a combination of b and a for an arbi- 
trary diagonal state f of & and a in the vacuum. Furthermore, 
since stochastic ordering is preserved under convex combina- 
tions, it suffices to prove the statement for any pure diagonal 
state f = fc 7^ 0. In this case the state of Catt is given 

by 



l=Q 



l + k 



\l + k){l + k\ 



1=0 



d\''>\l + k){l + k\ 



where we have relabelled ai by d and introduced the mode 
b = Rb^ + Ta\. As before, the state of b is the vacuum if 
and only if b is in the vacuum state, so it suffices to verify the 
statement for the state f = \k){k\ in which case the output 
state is 



-2g{k+ 



l + k 



10(^1 = Erff^io(/i 



The relation < now follows from 

p 



1=0 1=0 ^ ^ 



< 1 



1=0 

p 

1=0 



(0) 



This ends the proof of Lemma A. 2 for both cases 



□ 



The following lemma completes the proof of Theorem 11. 1 
by transforming the stochastic ordering into the desired norm 
inequality ^Al) . Its proof |[l4l|2Tl uses the fact that q ^ p'^ 
which is equivalent to the fact that P*(<I>*i) is more noisy 
that (f>*2 . The latter is satisfied for /c > fcp as assumed in the 
theorem. 

Lemma A.3. Let p' be a discrete probability distribution such 
that < p' . Then 



lb'-'?lli>lb^ 



(A4) 

□ 



where F — tanh(i) and — cosh(?). The stochastic order- 
ing now reduces to showing that Ym^q ^[^^ — J2]^o foi" 
all m. With the notation 7 = F^, we get 



p+k 



(=0 



k 



1 — n V / 



(0) 



1=0 



Amplification. As before we write a = cosh(t)ai + 
sinh(t)a2 and define i by cosh(?) = k cosh(i) and the beam- 
splitter coeficients 



sinh (t) 



sinh(t) 



Proof of LemmaliL2l 

We use the notations introduced in the proof of Theo- 
rem II. 1 By expressing the quadrature variance of the input 
mode a in terms of t and si we obtain sinh^ t = -^r-^- Ac- 
cording to Theorem |II.l| the output state of the optimal channel 
P* is the Gaussian state 



-2gsl 



y\i){i\ 



1=0 



with g taking different values in the attenuation and amplifica- 
tion cases. For the geometric distributions and q we have 



|$^-$^-^||i = ||p"-g||i = 2(s 



mo + 1 



where mo is the largest integer such that p':^^^ < q„ig, more 
precisely 

mo= Lln[(l~S)/(l-S2)]/ln(.S2/S)J 
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It remains to compute the concrete expressions of .1 and im- 
plicitly of mo for the attenuation and amplification cases. For 
attenuation, making use of sinh i — k sinh t, we find 



1 — Si + Si/C^ 

For amplification, we use cosh t — k cosh t and find 



^amp 



= 1 



A;2 ■ 



Proof of Theorem lllLTl 

Proof. We want to show that Q* := Sn° P* o Tn is the opti- 
mal purification or dilution procedure for n i.i.d. qubits. The 
idea is that, by using LAN, we can show the qubit and Gaus- 
sian statistical problems to be equivalent, the Gaussian one 
(respectively for attenuation and amplification) being solved 



in Section 11 B which then allows us to recast the qubit prob- 



lem in the Gaussian setup with a vanishing difference in the 
risks. We will consider the four separate cases: zero risk, 
solely quantum contribution, solely classical contribution, and 
both classical and quantum contributions. We will then use the 
Gaussian solution to show that i?max(roi Q*i ^) is less than or 
equal to the corresponding optimal Gaussian risk, then show 
that a strict inequality violates the optimality of this optimal 
solution. We begin by restricting r to the local neighbourhood 
I |r — Toll < n^2+^. This probability that the state fails to be 
in this region is o(l) and has no influence on the asymptotic 
risk (see Lemma 2.1. in ifTOl ). We are now able to apply LAN, 
which maps input states pf" close to some Gaussian state, say 
$„, via the channel T„. We now consider the individual cases, 
which are each slight variations on the same proof: 

Casel: k < min(fco, fcp'^'*)- In this case both classical and 
quantum Gaussian channels have zero risk, so the asymptotic 
qubit risk is zero. 

(c) 

Case 2: fcp < fc < fcp . In this instance, the risk receives 
only a quantum contribution. Using contractivity of the CP 
maps Sm and P*, the LAN convergence, and the fact that in 
this regime K*Nu ~ N^' we obtain 

R{Qn,r,X) 

= \\p^-Qn{Pu)h 

+ \\SmiK' ® ^u')) - SmiP-" ® K*{T„{p:)))\\l 

< \\K' ^ ^u' ~P*(^ K^{TM)\\i + o{l) 
= \\K' ® ^u' -P*® i^*($u' ^ ^u)||l + o(l) 
= ||4>^?-P*($^^)||i+o(l) 

= -Rmmmax(si,S2,fc) +o(l) (A5) 

where i?ininmax(si, S2, k) is the minimax risk for the quan- 
tum Gaussian problem, obtained in Theorem II. 1 By taking 
supremum over ||u|| < n*^ we get 



which implies that 

-^minmax ('"o 1 — -^minmax ("^1 ; ^2 5 

Next, we show by contradiction that this inequaUty cannot 
be strict. Suppose that there exists a sequence of purification 
or dilution procedures Qn, which act on qubits and satisfies 

^max((9nro, A) < i?ininmax (si , S2 , fc)"?? for SOmC 7/ > and 

n > uq. We will use LAN to show that there exists a Gaussian 
dilution (amplification) channel whose risk is strictly smaller 
than the minimax risk, which is a contradiction. 
The general setup can be seen in (|A6|l 



^ 11 



(A6) 



Here LAN is restricted to a two dimensional family of rotated 
qubit states, for which the limit model is quantum Gaussian, 
with no classical component. Assuming is in the domain 
of applicability of LAN (which can be effected by an adaptive 
measurement lfT4l '). we get the inequalities 

\\T^oQ„oSn{K')-Kn\ 

<\\Qn{p:)-p'^\\i + \\TM)-Kn\i 

(si, S2, fc) - ?7 + 0(1). (A7) 



< R 



inmmax \ 



mi 11 max 



(si,S2,fc). 



By taking the limit n — > oo we get the desired contradiction. 

Case 3: /cq^^ < k < k^. This case applies only to dilution 
and the risk receives only a classical contribution. The proof 
follows the same steps as the previous case, with the quantum 
Gaussian replaced by the classical one. The inequality ( |A5[ ) 
becomes 

-IIpu'-Q:.(pu)IIi 

< ||7V(0, 1 - A^llroir) - X*(iV(0, 1 - llrof ))|li + o(l) 

- i?minmax(^l, ■^^2, fc) + o(l). (A8) 

where i?minmax(Vi, V2, fc) is the optimal risk of Eq. (BJ and 
we have identified Vi = 1 ~ ||ro|p, ¥2 = !- A^HroHrThis 
implies 

-^minmax (1*0 7 A) ^ -^minmax (^l,^2,fc). 

The equality is obtained by showing that strict inequality 
would lead to a classical amplification procedure whose risk 
is smaller than the minimax risk. 

Case 4: k > max(A:Q'^\ fco)- This case applies to both di- 
lution and amplification, and both the quantum and classical 
channels contribute to the risk. 

= IIpu"-Q:(pu)IIi 

<\\K''^ ^u' -K*®P* ($^1 ® iVu) II 1 + 0(1) 

< -Rminmax(14,V2,Si,S2,fc) +0(1). (A9) 
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Here i?minmax ( Vi , 1^2 , Si , S2 , fc) is the minimax norm-one risk 
for the problem of transforming the Gaussian state ^f,^ (E) 
into (E) N^'- Since the quantum and classical components 
are independent and have different local parameters, the opti- 
mal channel is the product K* P*. The explicit expression 
of the minimax risk is 



-Rminmax(Vi, 1/2, Si, S2, fc) 
f 



exp 



'27rfc2(i_||r„||2) 

2(l-A2||ro|P) 



(1 - 



2^(l-A2||ro| 



(1-S2)s^ 



(AlO) 



Finally, the equality in ( A9 1 can be proven by contradiction as 
in case 2. □ 
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